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Abstract: 

We present a supersymmetric model with two dark matter (DM) components explaining 
the galactic positron excess observed by PAMELA/HEAT and ATIC/PPB-BETS: One is 
the conventional (bino-like) lightest supersymmetric particle (LSP) Xi ^md the other is a 
TeV scale meta-stable neutral singlet No, which is a Dirac fermion [N, N^). In this model, 
No decays dominantly into xe^e~ through an R parity preserving dimension 6 operator 
with the life time tn ~ 10^^ sec. We introduce a pair of vector-like superheavy SU(2) lepton 
doublets {L,L^) and lepton singlets {E,E^). The dimension 6 operator leading to the 
decay is generated from the leptophilic Yukawa interactions hy W D Ne^E + LhdE^ + 
m^/2hL'^ with the dimensionless couplings of order unity, and the gauge interaction by 
£ D \/2g' e'^* e'^X + li-C- The superheavy masses of the vector-like leptons {Ml, Me ~ 10^^ 
GeV) are responsible for the longevity of N^. The low energy field spectrum in this model 
is just the MSSM fields and Nd. Even for the case that the portion of Nu is much smaller 
than that of x in the total DM density [0(10"^*^) < njv^/n^], the observed positron excess 
can be explained by adopting relatively lighter masses of the vector-like leptons (10^^ GeV 
< Ml,e ^ 10^^ GeV). The smallness of the electron mass is also explained. This model 
is easily embedded in the flipped SU(5) grand unification, which is a leptophilic unified 
theory. 
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1. Introduction 

Although the existence of dark matter (DM) is advocated through several cosmological 
observations, we don't know yet its identity. The lightest supersymmetric particle (LSP) 
in the minimal supersymmetric standard model (MSSM) has been believed to be one of 
the most promising DM candidates over last two decades. Not only it is well-motivated 
from the promising particle physics model i.e. the MSSM, but also it naturally carries the 
features required for a weakly interacting massive particle (WIMP). 

However, recently the PAMELA group reported very challenging observational results 
on excess of high energy positrons coming from the galactic halo It has confirmed the 
previous similar observations of HEAT with the improved precisions. The PAMELA 
data shows a rising positron flux e^/(e^ + e~) from 10 to 100 GeV, which gives rise to 
considerably big deviations (~ 0.1) from the theoretically expected values Q. On the other 
hand, any significant anti-proton excess was not observed. Indeed, these observations would 
be quite hard to explain in large classes of models only with one DM component of the 
Majorana fermion such as the MSSM. Moreover, the ATIC and PPB-BETS collaborations 
also reported recently their observations that the flux of {e~^ + e~) keeps rising upto around 
800 GeV @, It might imply that if such deviations result from the DM's annihilation 
or decay, the mass of the DM would be in a TeV range. 

Various scenarios beyond the conventional DM models have been suggested so far, 
including new scenarios of DM annihilation [|, |, |], DM decay |, |l|, 0], and nearby 



pulsars producing charged leptons [|l^]. However, if the mass of DM is above TeV, and if 
we keep the galactic profile of NEW and Einasto [l^, the annihilation scenario would be 
disfavored due to the bound of the 7 ray flux by the HESS observations of the galactic 
ridge [^]. (See also the recent discussions on DM annihilation in Refs. p^].) 

In Refs. 1^ 0], the PAMELA/HEAT anomaly could be explained by co-annihilations 
of the LSP and another DM component No, which is a Dirac fermion (N, N^) with a weak 
scale mass. It turns out that in this case, introductions of a pair of extra vector-like lepton 
singlets {E, E^), and a leptophilic coupling of the DM, Ne^E are indispensable to explain 
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the PAMELA/HEAT's observations. The relatively small masses of {N,N'') and {E,E'') 
[~ 0(100) GeV] and the needed coupling in the superpotential for decay of E and 
into the DM can be guaranteed by the U(l)i? symmetry. This model ('WdmMSSM") is 
easily embedded in the flipped SU(5) grand unified theory (GUT), which is assumed to be 
broken to the MSSM at the GUT scale 0, 0]. 

In this paper, we attempt to explain the anomalies from PAMELA/HEAT and ATIC/PPB- 
BETS by decay of the extra DM component No- If the observations by PAMELA/HEAT 
and ATIC/PPB-BETS are caused indeed by DM decay, the data of ATIC/PPB-BETS 
imply that the mass of the DM matter would be around 2 TeV and its life time ~ 10^^ sec. 
Such a long life time is achievable, if the TeV scale DM decays to (+ neutral particles) 
through a dimension 6 operator suppressed by the mass squared of order the GUT scale 
|1C, 11]. To be consistent with the PAMELA data, the hadronic decay modes of DM should 
not exceed 10 % |0|. 

This paper is organized as follows: In Sec. 2, we discuss how a dimension 6 operator 
for DM decay can be naturally dominant in a supersymmetric model extending the MSSM, 
and in Sec. 3, we propose a specific model realizing the idea discussed in Sec. 2. In Sec. 4, 
we briefly discuss the implications of the recently observed anomalies in phenomenological 
model building in string theory, and summarize our conclusions. 



2. Dark Matter Decay 

For the case that DM decays to the positrons, the positron flux at the earth is given by a 
semi-analytic form [p!c|]: 

$e+(i?) = ^ Tdm X -— / dE' (2.1) 
\mj)Mj 4h[E) Je dE' 

where dN^+fdE' is the spectrum of e+ produced by DM decay, p the DM energy den- 
sity at the earth, and mDM its mass. b{E) [= £'^/(GeV-10^^ sec.)] indicates the energy 
loss coefficient, and I{Xd) is the halo function for DM decay, which depends only on the 
galactic astrophysics. As mentioned in Introduction, if DM decay causes the observations 
by ATIC/PPB-BETS, it implies the DM mass moM should be about 2 TeV. For p w 0.3 
GeVcm"^, the decay rate Fdm should be around 10~^^ sec.~^ This value can be achieved if 
the DM decays to light leptons through a dimension 6 operator mediated by a GUT scale 
massive field (~ 10^^ GeV) ]10|: 

It might imply that the observations of ATIC/PPB-BETS could be interpreted as signals 
of the GUT scale physics. 

If DM with a TeV scale mass was the LSP such as the neutralino or gravitino, and so 
supersymmetry (SUSY) breaking soft masses of the visible sector fields should be heavier 
than TeV, then the status of SUSY as a solution of the gauge hierarchy problem in particle 
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physics becomes weak and a considerable fine-tuning for the Higgs mass would be unavoid- 
able. Thus, we introduce a Dirac type extra DM component {N,N^} (= Nd) with a TeV 
scale mass apart from the (bino-like) LSP, x- If ^ ^■iid N'^ are the dominant component 
of DM, the mass of the heavy field mediating DM decay should be around the GUT scale. 
However, we have one more DM component Xi which is the absolutely stable particle. If 
the portion of x in the total DM density is large (and so the portion of and is small), 
the mass of the heavy field mediating DM decay needs to be adjusted such that the fiux 
needed to explain the positron excess is fixed. That is to say, in Eq. (|2.1| ) the smaller DM 
number density by a factor (M^/Mg)^ (< 1) can be compensated with a larger decay rate 



by assuming the lighter mediator in Eq. ( ^.21) such that the Eq. {2A) remains the same: 

M \^ 

for Mg^M, in Tdm- (2.3) 



TTlBM "IDM \Mg 

However, we should protect the results of the standard big bang necleosynthesis. For the 
life time of N longer than the age of the universe I/Fdm > 10^^ sec, the ratio of the 
number density of to X) ^A^/^-x is indeed extremely flexible: 

O(lo-iO) < ^ for 10^3 GeV < Af, < 10^^ GeV. (2.4) 

Even with an extremely small number density of N^, thus, the energetic positron excess 
from the recent experiments can be easily explained. 

In fact, how much Nd was created in the early universe is quite model dependent. In 
Ref. 1 11], when DM interacts with the standard model particles via a GUT suppressed di- 



mension 6 operator, the reheating temperature to produce DM much enough to explain the 
energy density pdm ~ 10~^ GeVcm"^ is estimated as 10^° GeV. Even if the reheating tem- 
perature is much lower than 10^'^ GeV, there exist many other possibilities to produce No 
sufficiently, depending on inflationary scenarios. For instance, could be non-thermally 
created directly from the inflaton decay or indirectly via a hidden sector fleld decay. Or 
A'^ could be coupled to the hidden sector flelds X and X"^ with TeV scale masses through 
W D NXX^. Then A^, A^'^ could be in a thermal equilibrium state with X, X (and also 
X^, X'^) by exchanging their scalar partners X'^ {X) down to a proper decoupling temper- 
ature, which is deflned with hidden sector flelds. However, we do not specify a possibility 
in this paper, because we have two dark matter components and so have extremely large 
flexibility for the portions of UN/n^, as mentioned above. 

As in the co-annihilation DM scenario of Refs. [^, 0], the electrophilic coupling of DM, 

Ne'^E (2.5) 

is essential also in the DM decay scenario to be consistent with the PAMELA'S observations: 
They did not observe the anti-proton excess from the cosmic ray. Here indicates the flrst 
family of the lepton singlet in the MSSM, and E is a newly introduced heavy lepton singlet 
with the same electric charge as e~ . To cancel the anomaly, E should be accompanied with 
E'^, whose electric charge is opposite to that of E. They achieve heavy masses (of order 
the GUT scale) from the Dirac mass term MeEE^. 
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Of course, one could think the possibihty that the singlet couples also to the quark 
singlets and extra vector-like heavy quarks. In this case, however, the coupling with them 
should be relatively small such that the hadronic decay rates do not exceed 10 % for the 
consistency with the PAMELA'S observations p!o| ]. In this paper we do not introduce such 
heavy extra quarks for simplicity of our discussion. 

Since the A^ is required to discriminate leptons and quarks, the interaction Ne'^E can 
not be accommodated in the conventional GUT such as SU(5) and SO(IO). It is, however, 
well embedded in the flipped SU(5) [= SU(5)xU(l)x] 0, which is a phenomenologically 
promising GUT. Indeed, flipped SU(5) is the leptophilic GUT. It means that the leptons, 
particularly the lepton singlets are special in flipped SU(5). Since the lepton singlet e'^ 
remains an SU(5) singlet (= I5) under SU(5)xU(l)x, the A^ (= Iq) can couple to and 
E (= I-5) without being accompanied with quarks, and Ne'^E (= I0I5I-5) is invariant 
under the flipped SU(5) gauge symmetry. 

For dimension 6 dominance in the A^ decay, the dimension 4 and 5 operators from 
Nhuhfi, Nlihu, Nlihde'^/Mp, Nlilje'^/Mp, etc. in the superpotential, where hu, h^, and k 
{i = 1,2,3) indicate the MSSM Higgs and lepton doublets, should be removed from the 
Lagrangian by introducing a proper symmetry, because they open too fast A^'s decay chan- 
nel into X s-nd the standard model leptons (r^v ~ m^/Mp ~ 10~^ sec.~^, if the dimension 

5 operators are dominant.). The dominant GUT suppressed dimension 6 operator would 
be generated from renormalizable operators, in which GUT scale heavy fields are involved. 
The Feynman diagram displaying such induced dimension 6 operators would satisfy 

(Number of heavy fermion propagators) — (Number of heavy mass insertions) = 2. 

(2.6) 

We introduce a pair of vector-like heavy SU(2) lepton doublets (L, L'^) with their Dirac 
mass term M^LL'^, where ~ Mq ~ 10^^ GeV, in the superpotential for the dimension 

6 process. 

In SUSY theories, a diagram replacing some heavy fermions lines by the scalar partners' 
lines is always present. The heavy scalars' mass squareds are the same as the fermions' 
upto the SUSY breaking soft mass squareds of ©(rrigyg)- Since the contributions by the 
scalar propagators are suppressed by ~ 1/Mq at low energies, the diagram by the heavy 
fermions dominates over the diagram replacing some of them by their scalar partners. 

If the couplings between E^ and the MSSM lepton doublets li such as lih^E'^ (and 
also liljE^) were present in the superpotential, Nlihde^ /Me (and Nlilje^ /Me) could be 
induced after integrating out E and from lih^E'^ (liljE^) and Ne'^E. Therefore, lih^E^ 
and liljE^ also should be disallowed from the superpotential by a symmetry. Even though 
they were generated when the symmetry, which forbids them, is broken, they are still safe 
only if their coefficients are suppressed by 0(111^/2/^^0) ■ 

On the other hand, we require the presence of 

LhdE" and my2hL'' (2.7) 

in the superpotential in order to make it possible for A^ (and also A^'^) to decay into x^^e^ 
via the dimension 6 operator. Here h stands for the first generation of the lepton doublet 
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Figure 1: Dominant diagram of N (and iV^) decay: The dimensionless Yukawa couplings are of 
order unity. The mass parameters Me and Ml are of order the GUT scale. 



in the MSSM, and indicates a TeV scale mass parameter. See the Feynman diagram 
in Figure 1 for the process 



X + e + 



(2.8) 



It is the dominant diagram of N and decay, if mjv < . Its decay rate is estimated 
as 



m 



N 



N 



1927r3 



2m%MEML 



(2. 



where Fat ~ IQ-^^ sec.-^ for ~ 2 TeV [> 10 x 0{m^)], Me ^ Ml ^ 10^^ GeV, and 
the contributions by the dimensionless Yukawa couplings 0(y^) ~ 1. 

If the selectron e*^ is relatively light, m^v ^ mg':, however, the decay channel N — > 
+ e^* is kinematically allowed, and e^* can be an on-shell particle in Figure 1. Then, 
the decay rate becomes enhanced by 0(100): 



[m 



N 



N 



m 



2 \2 



j'{hd)m3/2 
MeMl 



(2.10) 



For Fjv giving 10^^^ sec^"*^, however. Me ~ Ml ~ 10"*^^ GeV is not much affected. 

The N (and also N'^) could decay also into the hadrons through the Higgs of Figure 1. 
However, such decay modes are much more suppressed than 10 %. It is because they are 5 
body decay channels with much suppressed phase factors, and the quarks' (except the top 
quark) Yukawa couplings at the vertices, from which quark branches start, are so small. 

L, and are embedded, respectively, in 5_3, 63 and 5_2 in flipped SU(5), which 
are accompanied with quark singlets. Since Ne'^E (= I0I5I-5) is still electrophilic and 
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the doublets/triplets in 52, 5_2 are split by the missing partner mechanism |17|, however, 



the decay channels to hadrons are suppressed also in the flipped SU(5) model. ^ 
3. The Model 

The relevant superpotential is composed oi W = VF/vdecay + W^mass, where Wj^decay and 
^mass are, respectively, given by 

WNdec^y = Ne'^E + LhaE"^ + + mg/s/iL^ (3.1) 
H^mass = MlLL^ + MeEE' + TUnNN', (3.2) 

where we drop the dimensionless Yukawa coupling constants for simplicity. All the dimen- 
sionless Yukawa couplings in the above expressions are tacitly assumed to be of order unity. 
Not introducing N couplings to quarks, we can avoid anti-proton excess. We assume that 
the masses of L, L^, and E, E"", i.e. Ml and Me are of order Mq (~ 10^*^ GeV). On 
the other hand, the masses of the DM N, N'^ should be constrained to be about 2 TeV in 
order to explain the observations of ATIC/PPB-BETS. In fact, the mixing term liL^ and 
the mass terms of N, and also h^, all of which are proportional to ^1,3/2, can not 
be present in the bare superpotential due to the U(l)ij symmetry, which will be presented 
later. However, via the Giudice-Masiero mechanism |jl8|, the Kahler potential 

^(A/iL-^ + AXM+h.c. (3.3) 

can induce such supersymmetric mixing term seen in Eq. ( |3.1| ), and also the MSSM "^u 
term" , if SUSY is broken by a nonvanishing VEV of the F component of the singlet super- 
field S. We assume that (F^) ~ m3/2Afp, where m3/2 ~ TeV, A ~ 0{1), and A' ^ 0(0. 1). 
So we have the term of O(IOO) GeV. We will explain later how the Dirac mass term of N 
and A'^'^ in Eq. ( |3.2| ) (ttzat ~ 2 TeV) is naturally generated. Concerning the gaugino mass 
terms, the gauge kinetic functions are involved in their mass generations in supergravity. 
We assume that the LSP is the (bino-like) neutralino with the mass of O(IOO) GeV. It is 
stable and also a component of the DM together with Nd. 

The cubic term of N in Eq. (|3.l| ) is introduced such that the superpartner of N, i.e. 
N promptly decays to 2 A. We just assume that the soft mass of A^ is heavy enough (> 4 
TeV) for this decay process to open. Were it not for A^^ in the superpotential, -A could 
remain meta-stable together with A^. Actually it is not a serious problem. But we prefer 
smaller number of species of DM. By the SUSY breaking B-term corresponding to the third 
term in Eq. ( |3.2| ), A^ can be converted to A, and so jSf^ also can decay to 2 A through the 
A'^ term. 

^In flipped SU(5), lOi and 5-3 contain {cr,q,i>'^} and {u'^,1}, respectively, where d'^ {u'^) denotes the 
quark singlet of Qcm. = 1/3 (—2/3), and q (l) is the quark (lepton) doublet. Note that the SU(2) singlets 
are interchanged by each other, compared to the Georgi-Glashow's SU(5). Particularly, while the Majorana 
neutrino v'^ is included in the tensor multiplet, the charged lepton singlet e'^ {— I5) is the SU(5) singlet. 
The MSSM Higgs doublets hu and ha are embedded in 62 (= 5^) and 5-2 (= 5h), respectively, in flipped 
SU(5). The flipped SU(5) gauge group is broken to the standard model gauge group by the Higgs flelds of 
the tensor representations IOh, IOh- In flipped SU(5), the doublets/triplets included in 5h, 5h are simply 
split just through lOnlOfrSh and lOfrlOnSh. 
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Table 1: The quantum numbers of the superfields. The superfields written with the capital letters 
are the newly introduced fields, which are absent in the MSSM. Except A'', N'^ and the MSSM 
Higgs fields, the vector-like superfields are all decoupled from low energy physics due to their heavy 
masses. 



The global symmetry observed in this model is U(l)ijxU(l)pQ x Z™, where de- 
notes the matter parity (or R parity). The quantum numbers of the superfields under the 
symmetry are displayed in Table |l|. It is straightforward to assign the quantum numbers 
also to all other MSSM chiral superfields, which are not presented in Table so as to 
admit all the needed R parity preserving Yukawa terms in the superpotential. 

violate the \J{1)r and U(1)pq symmetries. 



In fact, the last terms of Eqs. (3.1) and (|37 
They and soft SUSY breaking A- and B-terms corresponding to the superpotential of 
Eqs. (3d.) and (|3.2| ) in the scalar potential break the R symmetry to Zq. That is to say, 
SUSY breaking effects result in U(l)/j breaking into Zg. Since S carries the unit charge of 
U(l)pQ symmetry, it is also broken at the intermediate scale (~ y^m^J^Mp ~ lO^'' GeV). 
A proper linear combination of the symmetries could still remain unbroken, but we will see 
that it is also broken by VEVs of other fields. 

Among the unwanted terms in the superpotential, which were discussed in Sec. 2, 
Nlihde'^ and lihdE'^ are induced from the bare Kahler potential: K D Ti^ {Nlhde'^ /Mp + 
IhdE^ / Mp) -\- h.c, because they are consistent with the charge assignments in Table |l|. 
But their suppression factors m^/2/Mp and m^i2/Mp are small enough. However, Nlihdc'^ 
and lihdE'^ suppressed, respectively, by m^/2/MQ and 711^/2/^0 rather than m^/2/Mp and 
171^12/ Mp can be generated. They are exactly what are shown from Figure 1. 

Indeed, the diagram in Figure 1 is reminiscent of that explaining the seesaw mecha- 
nism of the neutrino. Integrating out the superheavy fermions E, E^^ and L, L'^ yields 
the effective Lagrangian relevant for N — > e*^ + The equations of motion, dC/dE = 
dC/dE'^ = dC/dL = dC/dL^ = 0, where E^^\L^'^^ are the fermionic components of the 
corresponding superfields, give E'^ = —c^N/Me, E = —{hd)L/ME, L'^ = —{hd)E'^ /Ml, 
and L = —m^i2li/ML, respectively. By inserting them back into the original Lagrangian, 
one can get the effective Lagrangian, and also the effective Kahler potential: 



m 



eff. 



3/2 



MeMl 



hde%N C 



St 



MpMeMl 



hde^liN + h.c. 



(3.4) 



Thus, A'' (and N'^) can decay to e^e'^* with the suppressed amplitude ~ m^/2{^d) /M^ 



2_ 



By the gauge interaction C D \pig'e*e^x + h.c. 
Figure 1. 



e*^* eventually decays to e x ^ shown in 
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Table 2: The quantum numbers of some superfields. They all are decoupled from low energy 
physics due to their heavy masses. 



(5) 

I 
I 

T 



Ml' 
— ^ 



(hd) 
I 

t 



h L"' L' 

Figure 2: Small electron mass generation: The dimensionless Yukawa couplings are of order unity. 



The charge assignments in Table |T] forbid the Yukawa coupling for the electron mass, 
lihde^ from the bare superpotential. After \J{1)r and U(1)pq broken, however, one can 
expect the electron mass term is generated. Let us consider the following superpotential: 

W^eiec = ShL^' + L'hde' + Ml'L'L^' + -^S^S^ + (3.5) 

ivip IVlp 

where the quantum numbers of S, S, and L', L'^ are presented in Table ^. The scalar 
potential derived from the last term of Weieo the A-term corresponding to it, and soft mass 
terms of S, S could allow the minimum, where nonzero VEVs of S and S are developed 
(but (N) = {N'^) = 0), breaking the U(1)r and U(1)pq symmetries completely: 



(S) ~ (S) ~ ^m^/^Mp ~ 10^" GeV. (3.6) 

Namely, SUSY breaking triggers the PQ symmetry breaking at the same energy scale. 
However the symmetry is still unbroken, which plays exactly the role of the matter 
parity (or R parity) in the MSSM. The Zf^ parity conservation still prevents the dangerous 
term liljE'^ as well as liljC^ from being induced. Even with the nonzero {S) and (5), the 
other unwanted terms also don't appear at low dimensions. 

Once S develops a VEV, the desired mass term of DM N and N'^ can be achieved via 
the last term of Eq. (|3^ ): 

mN = ~ 2 TeV. (3.7) 

Mp ^ ' 

The electron mass term {{S) /ML')lihde'^ also can be generated as seen in Figure 2. It 
dominates over the nonrenormalizable term in the bare superpotential, {{S) /Mp)lihde'^. 
Thus, the correct order of magnitude of the electron mass can be achieved if M/^/ ~ 10^^ 
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GeV: 



nip. 



Ml' 



10-6 



(3.8) 



Expanding S in terms of the axion field a fT|,|2^, S" = {FaNDW + p)e'"^'^^''^°'^\ where 
Fa (= {S)/N£)\Y ~ 10^'' GeV/NDw) is the axion decay constant and indicates the 

domain wall number, (S/Mx^f)lihiie^ yields the axion-electron coupling: 



FaNow 



(3.9) 



In fact, if there is a energy loss mechanism by a very weakly interacting light particle 
like the axion, it can affect the resulting luminosity function of the white dwarf. Thus, 
provided that the axion is dominantly involved in the cooling mechanism of white dwarfs, 
the luminosity function of the white dwarf can be used to estimate the axion decay constant 
Fa'. If the axions lighter than 1 eV are produced inside white dwarfs, the axion's coupling 
to the electron is estimated as 0.7 x 10~^^ for the best fit of lH) 0' 



meT{e) 



FaNow 



0.7 X 10 



-13 



where r(e) denotes the PQ charge of e. Thus, Eq. (|3.10| ) means FaN^w 
GeV for r(e) = ±1, which can be coincident with Eq. (^). 



(3.10) 



0.72 X 10^° 



4. Discussion 



The weak scale SUSY is important, because it provides the resolution of the gauge hierarchy 
problem and a way to connect the standard model at low energies and the string theory at 
high energies. If the low energy SUSY should be accepted as a basic language describing our 
nature, the DM in the universe should be understood within this framework. However, the 
recently reported observations on excess of energetic positrons from cosmic ray might be 
quite embarrassing, because it is seemingly hard to understand in terms of the conventional 
MSSM. 

In this paper, we present a SUSY model with one more dark matter component with 
a TeV scale mass apart from the LSP. The observed anomaly could be naturally explained 
by the extra DM's decay through a dimension 6 operator, which is induced by the renor- 
malizable operators. The life time of DM (~ 10^^ sec.) needed to explain the observed 
positron flux is caused by the heavy masses (~ 10^^ GeV) of the mediators involved in 
the decay process. Since this model permits the possibility that the field spectrum below 
the GUT scale can coincide with that of the MSSM except for No, the gauge coupling 
unification in the MSSM could be still maintained. Therefore, we could rescue the string 
models realizing the MSSM with siv?Ow = 3/8 ||2^, because a lot of neutral singlets are 
easily found in string models. 



-9- 



Acknowledgment s 

The author is supported by the FPRD of the BK21 program, in part by the Korea 
Research Foundation, Grant No. KRF-2005-084-C00001 and the KICOS Grant No. 
K20732000011-07A0700-01110 of the Ministry of Education and Science of Repubhc of 
Korea. 

References 

[1] O. Adriani et al. [PAMELA Collaboration], Phys. Rev. Lett. 102 (2009) 051101 
[arXiv:0810.4994 [astro-ph]]; Nature 458 (2009) 607 [arXiv:0810.4995 [astro-ph]]. 

[2] S. W. Barwick et al. [HEAT Collaboration], Astrophys. J. 482 (1997) L191. 

[3] L V. Moskalenko and A. W. Strong, Astrophys. J. 493 (1998) 694 [arXiv:astro-ph/9710124]. 

[4] J. Chang et al. [ATIC Collaboration], Nature 456 (2008) 362. 

[5] S. Torii et al. [PPB-BETS Collaboration], arXiv:0809.0760 [astro-ph]. 

[6] J. H. Huh, J. E. Kim and B. Kyae, Phys. Rev. D 79 (2009) 063529 [arXiv:0809.2601 [hep-ph]]. 

[7] K. J. Bae, J. H. Huh, J. E. Kim, B. Kyae and R. D. VioUier, Nucl. Phys. B 817 (2009) 58 
[arXiv:0812.3511 [hep-ph]]. 

[8] M. Cirein, M. Kadastik, M. Raidal and A. Strumia, arXiv:0809.2409 [hcp-ph]; 

N. Arkani-Hamed, D. P. Finkbeiner, T. Slatyer and N. Weiner, arXiv:0810.0713 [hep-ph]; 
M. Fairbairn and J. Zupan, arXiv:0810.4147 [hep-ph]; A. E. Nelson and C. Spitzer, 
arXiv:0810.5167 [hep-ph]; Y. Nomura and J. Thaler, arXiv:0810.5397 [hep-ph]; R. Harnik and 
G. D. Kribs, arXiv:0810.5557 [hcp-ph]; Y. Bai and Z. Han, arXiv:081 1.0387 [hep-ph]; 
P. J. Fox and E. Poppitz, arXiv:081 1.0399 [hep-ph]; S. Baek and P. Ko, arXiv:081 1.1646 
[hep-ph]; I. Cholis, G. Dobler, D. P. Finkbeiner, L. Goodenough and N. Weiner, 
arXiv:0811.3641 [astro-ph]; K. M. Zurek, arXiv:0811.4429 [hcp-ph]; E. J. Chun and 
J. C. Park, arXiv:0812.0308 [hcp-ph]; R. AUahvcrdi, B. Dutta, K. Richardson-McDaniel and 
Y. Santoso, arXiv:0812.2196 [hep-ph]; D. Hooper, A. Stebbins and K. M. Zurek, 
arXiv:0812.3202 [hep-ph]; X. J. Bi, P. H. Gu, T. Li and X. Zhang, arXiv:0901.0176 [hep-ph]; 
L Gogoladzc, R. Khalid, Q. Shafi and H. Yuksel, arXiv:0901.0923 [hep-ph]; W. de Boer, 
arXiv:0901.2941 [hcp-ph]; J. Hisano, M. Kawasaki, K. Kohri, T. Moroi and K. Nakayama, 
arXiv:0901.3582 [hcp-ph]; F. Chen, J. M. Chnc and A. R. Frcy arXiv:G9G1.4327 [hcp-ph]. 

[9] P. f. Yin, Q. Yuan, J. Liu, J. Zhang, X. j. Bi and S. h. Zhu, arXiv:0811.0176 [hep-ph]; 
K. Ishiwata, S. Matsumoto and T. Moroi, arXiv:0811.0250 [hep-ph]; arXiv:081 1.4492 
[astro-ph]; A. Ibarra and D. Tran, arXiv:0811.1555 [hep-ph]; C. R. Chen, M. M. Nojiri, 
F. Takahashi and T. T. Yanagida, arXiv:0811.3357 [astro-ph]; X. Chen, arXiv:0902.0008 
[hep-ph]. 

[10] E. Nardi, F. Sannino and A. Strumia, JCAP 0901 (2009) 043 [arXiv:0811.4153 [hep-ph]]; 
A. Arvanitaki, S. Dimopoulos, S. Dubovsky, P. W. Graham, R. Harnik and S. Rajendran, 
arXiv:0812.2075 [hep-ph]. See also K. J. Bae and B. Kyae, JHEP 0905 (2009) 102 

[arXiv:0902.3578 [hep-ph]]. 

[11] K. Hamaguchi, S. Shirai and T. T. Yanagida, arXiv:0812.2374 [hcp-ph]. 

[12] D. Hooper, P. Blasi and P. D. Scrpico, JCAP 0901 (2009) 025; H. Yuksel, M. D. Kistler and 
T. Stanev, arXiv:0810.2784 [astro-ph]; S. Profumo, arXiv:0812.4457 [astro-ph]. 



-10- 



[13] J. F. Navarro, C. S. Frcnk and S. D. M. White, Astrophys. J. 462 (1996) 563; 

A. W. Graham, D. Merritt, B. Moore, J. Diemand and B. Terzic, Astron. J. 132 (2006) 2685; 
J. F. Navarro et al, arXiv:0810.1522 [astro-ph]. 

[14] G. Bertone, M. Cirelli, A. Strumia and M. Taoso, arXiv:0811.3744 [astro-ph]. 

[15] F. Aharonian et al. [H.E.S.S. Collaboration], Phys. Rev. Lett. 101 (2008) 261104 
[arXiv:081 1.3894 [astro-ph]]. 

[16] P. Meade, M. Papucci and T. Volansky, arXiv:0901.2925 [hep-ph]; J. Mardon, Y. Nomura, 
D. Stolarski and J. Thaler, arXiv:0901.2926 [hep-ph]. 

[17] S. M. Barr, Phys. Lett. B112 (1982) 219; J.-P. Derendingcr, J. E. Kim and D. V. 

Nanopoulos, Phys. Lett. B139 (1984) 170; I. Antoniadis, J. R. Ellis, J. S. Hagelin and 

D. V. Nanopoulos, Phys. Lett. B 194 (1987) 231; B. Kyae and Q. Shafi, Phys. Lett. B 635 
(2006) 247. 

[18] G. Giudice and A. Masiero, Phys. Lett. B206 (1988) 480. 

[19] J. E. Kim, Phys. Rev. Lett. 43 (1979) 103; M. A. Shifman, V. I. Vainstein, V. I. Zakharov, 
Nucl. Phys. B166 (1980) 4933. 

[20] M. Dine, W. Fischlcr and M. Srednicki, Phys. Lett. B104 (1981) 199; A. P. Zhitnitskii, Sov. 
J. Nucl. Phys. 31 (1980) 260. 

[21] J. Isern, E. Garcia-Berro, S. Torres, and S. Catalan, Astrophys. J. 682 (2008) L109. 

[22] See, for instance, J. E. Kim and B. Kyae, Nucl. Phys. B770 (2007) 47 [hcp-th/0608086]; J. 

E. Kim, J.-H. Kim and B. Kyae, .JHEP 0706 (2007) 034 [hep-ph/0702278]; J. E. Kim and 

B. Kyae, Phys. Rev. D 77 (2008) 106008 [arXiv:0712.1596 [hep-th]]. 



- 11 - 



